Origin of branch points in the spectrum of PT-symmetric periodic potentials 
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There exists multiple branch points in the energy spectrum for some PT-symmetric periodic po- 
tentials, where the real eigenvalues turn into complex ones. By studying the transmission amplitude 
for a localized complex potential, we elucidate the physical origin of the breakdown of perturbation 
method and Born approximation. Mostly importantly, we derive an analytic criteria to determine 
why, when and where the bifurcation will occur. 
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I. INTRODUCTION 



Non-Hcrmitian (NH) quantum mechanics (QM) has 
been used to describe an open system because complex 
eigenenergies imply the existence of a sink or source for 
the particle. The PT-symmetric hamiltonian is a special 
case in NHQM that can exhibit entirely real spectra[l|. 
Its development has gone through various stages Q with 
applications ranging from quantum cosmologyj3j , quan- 
tum field theoryQ to quantum computation [J [||. A 
flourishing example is in connection with the optical lat- 
tice which can be simulated by the Schrodinger equation. 
This provides a fertile ground to test the NH-related con- 
cepts experimentally. In this spirit scientists have studied 
the PT-symmetric periodic potentials and investigated 
several interesting features: (1) The existence of sev- 
eral branch points which demarcate the real and complex 
spectraQ- (2) Optical soliton solutions are found within 
the energy gap with real eigenvalues Q . (3) A quantum 
phase transition is identified when a sufficiently large real 
periodic potential is superimposed to render all the spec- 
trum realQ. The latter two, which provide the basis 
for the fruitful beam dynamics [1, Q, are closely related 
to the branch point and bifurcation feature. Although 
the branch points and a complex spectrum have been 
attributed to the failure of perturbation expansion [ic|. 
we believe it is worthwhile to study in more detail the 
physical origin of this failure. 

We shall start by studying a localized imaginary po- 
tential before embarking on the more interesting periodic 
case. Several new features are discussed in Section II 
which include (1) Existence of an upper bound for the 
imaginary part of eigenenergies. (2) Breakdown of the 
Born approximation for scattering states unless a large 
enough real potential is added. (3) Fano resonance within 
the imaginary barrier can cause the transmission ampli- 
tude to become huge by repetitive enhancements. (4) 
An incoming wave packet will also exhibit special fea- 
tures when its momentum is near the resonance. We find 
that the size of the transmission amplitude is the key to 
features (2) and (3). In Section III we argue that the 
spectrum of the PT-symmetric periodic potential is pre- 
dominantly shaped by the transmission amplitude for a 
single unit cell. Therefore, the conclusions in Section II 



for a localized potential are still applicable in the periodic 
case. By use of the Bloch's theorem, we derive an ana- 
lytic criteria for the occurrence of branch point and the 
band structure bifurcation. Discussions and conclusions 
are arranged in the final Section IV. 



II. SINGLE IMAGINARY SQUARE POTENTIAL 
A. Localized states 

At first glance, one might think that the eigenenergies 
of a complex potential are most likely also complex. This 
is in fact not true. For instance, take an imaginary square 
potential iVo from x = to x = a where Vq > and 
zero potential elsewhere. It exhibits both real/imaginary 
eigenvalues, which are continuous/discrete and corre- 
spond to scattering/localized states, respectively. 

Since the barrier potential is symmetric with respect to 
x = a/2, the spatial part of the localized wavefunction, 
ip(x,t) = (j>(x)e~~ lEt where H has been set to unity for 
convenience, can be expressed as: 

±e -(ifc-«)(x-o/2) ; x<0 7 
A ^p(x-a/2) ± e -ip{x-a/2)y Q < X < a, (1) 
e (ik-q)(x-a/2) ^ X > a 

where momentum k and q are real and positive, p is com- 
plex, and the ± sign denotes respectively the even/odd- 
parity state with respect to x = a/2. The minus sign in 
the exponent of <p(x < 0) is to prevent it from blowing up 
at x — ► — oo. These momenta are related by the complex 
energy E: 



E 



-(ik — q) 2 p 2 
~ ~ 2m 



2m 



iV . 



(2) 



One should be aware that the standard concept of den- 
sity current need not survive ll( the transition to NHQM. 
For instance, in the case of a PT-symmetric potential 
j(x,t) has to be redefined as — i[ip*(x,t)dip(— x,t)/dx — 
ip(— x, t)dip*(x, t)/dx\ /2m in order to obtain the conti- 
nuity equation. For the sake of arguments, let us stick to 
the standard notion and deduce the following mathemat- 
ical equation from the Schrodinger equation for a general 



2 



complex potential V(x): 







d 



2ImV(x)W(x,t)\ 2 = - \i>(x,t)\ +—j (x,t) (3) 



dx" 



where the symbol Im denotes taking the imagi- 
nary part and jo(x,t) = —i\"ip*(x,t)dip(x,t)/dx — 
ifj(x,t)dif}*(x,t)/dx~\/2m. Integrating over the whole 
space gives 



dx^{x,t)\ 2 = 2 



dxlmV(x)\ip(x,t)\ 2 (4) 



where jo(x, t) vanishes at x = ±oo for the localized state 
in Eq.|T|). This leads to the following inequality 



dx\ip(x,t)\ 2 < Max(ImF) / dx\ip(x,t)\ 



ImE 



(5) 

which tells us that the height of the imaginary barrier Vq 
serves as an upper bound for ImE: 



< ImE < Vr 



o- 



(0) 



By matching the boundary conditions, we can further 
show that the number of these localized states is finite 
and their energies are discrete. In contrast to being zero 
for a real well, the current density j(x,t) is not only 
nonvanishing, but also quantized for an imaginary bar- 
rier/well: 



jn (x,t) = ^e- 2 ^ 



-a/2|+2Im_E„t 



(7) 



where n labels these discrete localized states. 

All the derivations so far are based on the assump- 
tion that Vo > 0. However, it is easy to generalize 
to an imaginary well. The action of changing the po- 
tential from iVo to — iVo is equivalent to reversing the 
time, which transforms the hydrant-like bound state into 
a sink. The same conclusion applies to the eigenfunctions 
and eigenenergies after we take the complex conjugate of 
the time-independent Schrodingcr equation. As a result, 
the upper bound in Eq.® needs to be modified as 



-V < ImE < 0. 



(8) 



where subscript r is to distinguish this right-moving state 
from its degenerate state </>;,fc(x) that travels in the oppo- 
site direction. Fitting the boundary conditions, we find 



t = 



-Akk'e 



-ika 



e i(k+k')a^fc _ f.l\2 _ e i(k—k')a(h _|_ £') 2 

2i(k 2 - k' 2 )sm(k'a)e- tka 

e i(k+k>)a( k _ k ,yi _ e i(k-k>)a( k + fc /)2 



(10) 
(11) 



where k! = p — iq is the complex momentum inside the 
barrier. The amplitude t in Eq.@ consists of contribu- 
tions from all paths that reflect inside the barrier arbi- 
trary number of times: 



t = t\t<xe 



k a 



t 1 r 2 r 1 t 2 e 3ik ' a + t 1 r 2 r x r 2 r 1 t 2 e' lk ' a 



(12) 

where £i j 2/7*1, 2 denote the transmission/reflection ampli- 
tudes at the boundaries x = and x = a, respectively. 
The higher-order terms in Eq. p2p converge to zero for 
a real barrier. But this is no longer true for a posi- 
tive imaginary barrier because the longer they stay in 
the source the more probability they can gain. As a re- 
sult, the values of resonance peaks can be much bigger 
than unity, and the Born approximation that neglects 
the higher-order terms is expected to break down. These 
two features are clearly shown in Fig[T]and the number of 
resonances and their peak values increase when we raise 
Vq in FigH 




EM 



FIG. 1: Set the effective mass of electron as m = 0.5. The 
transmission and reflection coefficients of an imaginary barrier 
with height Vb = 5 and width a = 2 are plotted as a function 
of E/Vo where E represents the energy of the scattering state. 



B. Scattering states 

Another interesting problem is to study the transmis- 
sion and reflection coefficients for a complex square bar- 
rier. This is relevant to our discussion of the spectrum for 
a periodic potential in Section III. The scattering state 
takes the form: 

( e ikx +re -ikx^ x < Qj 

<l>r,k{x) = I Ae l{ P~^ x + Be-^P-^, < x < a, (9) 
[ te lkx , x > a 



Another important issue of the scattering process is the 
interference condition for Fano resonance. In real barri- 
ers, the resonance in transmission results from the con- 
structive interference between different scattering paths 
shown in Eq. (|12[) . which gives the familiar condition 
a = n\/2 = rni/p. In the mean time, reflection peaks 
are attributed to the destructive interference when 

a={n+\)\ (13) 

is satisfied. These two conditions become more compli- 
cated for an imaginary potential because the phase shift 
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FIG. 2: Same as Fig[T] except the barrier height is increased 
by eight folds to Vo = 40. 



E 


223.6 


257.4 


294.1 


333.4 


\t\ 2 


247.5 


2322.5 


24485.9 


1507.7 


n 


13.52 


14.49 


15.48 


16.47 



TABLE I: Partial list of the energies and peak values of the 
transmission coefficient in Fig[2] and their corresponding n 
defined by n — 2a/\. 



due to reflections from the boundaries between the imagi- 
nary and real potential regimes is energy-dependent. The 
amplitude n,2 in Eq. (fl"2"]) can be determined by studying 
the reflection by an imaginary step potential of height 
iV : 



k' 



ri,2 



1, E « V Q 

iVo/AE, £>V 



(14) 



which shows that there is an additional phase shift rang- 
ing from to 7r/2 as E increases. It is straightforward to 
write down an equation similar to Eq. (fl"2|) for the reflec- 
tion amplitude, r. Since both equations are geometric 
series with the same common ratio, t\T2 exp(2i/c'a), it is 
expected that the transmission and reflective resonances 
for an imaginary barrier should all follow the same second 
criterion, Eq . (fl"3"|) . at the limit E Vo- To verify these 
interference features, a partial list of the energies, peak 
values of the transmission coefficient and corresponding 
n in Fig[2]arc numerically calculated and summarized in 
Tab. HI For an imaginary well, Eci. (fT4"|) remains valid af- 
ter Vq is replaced by — Vo- However, the resonance peaks 
are no longer pronounced because the sink now draws in 
probability as the wave reflects in the well. 

In order to facilitate the discussion of a PT-symmetric 
periodic potential in Section III, let us put the square bar- 
rier and well side by side to form a simple PT-symmctric 
potential. The perturbation expansion can be shown to 
still break down as in the periodic caseQ. The simi- 
larity goes further that the perturbation in both poten- 
tials can be salvaged by incorporating a big-enough real 
potential. On the other hand, we can ask how strong 
a complex potential can ruin the perturbation expan- 



sion. Take this pair of square barrier and well potential, 
V(x) = {U [Q(x+a) -e(x-a)}+iV [e(x+a) -Q(x)]- 
iVo[Q(x) — Q(x — a))}, for example. The threshold for Vq 
is determined when the maximum of the common ratio 
in Eq. (fl"2"]) equals 1; namely, 



where 



n = 



'2 



|riT2 exp(2ik'a)\E= 



VE-VE-Up-iVo 
y/E + VE-Uo-iVo 
y/E-U -iVo-y/E-U + iV 
^E-U - Wo + y/E-U + Wo 



(15) 



k' = ^/2m{E -Uo-iVo). 



(16) 



For the parameters m = 0.5, a = 1 and Uo — 50, this 
critical Vo equals 0.00136. Numerical calculations show 
that this is the same threshold for the transmission am- 
plitude to exceed 1, which shall be proved in Section III 
to be the harbinger for the bifurcation of specturm for 
the PT-symmetric periodic potential. 



C. Evolution of a wave packet through the 
imaginary barrier 

When a wave packet impinges on a real potential bar- 
rier, the resonances are not critical to its evolution since 
their transmission amplitude is of the same order as the 
nonresonant ones. However, as is shown in Tab. HI the 
resonant amplitudes can reach a few tens of a thousand 
for an imaginary barrier. It is, therefore, of interest to 
study how a wave packet progresses through such a po- 
tential. 

As usual, we first decompose it into the scattering 
states 4>k in Eq.© and the localized states <p n in Eq.([T]): 



*(x,t = 0) = J2C k ct>K(x,0) 



E 



C n tf(x,0) (17) 



where the superscript R denotes the right eigenstate. 
This is similar to the case of nonHermitian matrices, for 
which the right and left eigenstates may not be the same. 
Since they can exhibit different physical properties and 
time evolution, it is important to distinguish them. Take 
the imaginary square potential, for instance. The barrier 
seen by the right eigenstates in Eqs.© and (JTJ) will turn 
into a well for their left counterparts, which are defined 
by 



\H = 



\E 



k/i 



or 



f P 
1 2m 



iV [@{x) - e(x - a)]\cf> 



77"* 



ft,* 



(18) 



(19) 



when expressed in the real space. These two eigenstates 
obey the following orthogonality relation 

(4>k/ n \4>k' /n') = (^fc/nl^fc/nMfc/n.fc'/n'- 



(20) 
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FIG. 3: The evolution of the wave packet in Eg. (|22[) impinges 
on an imaginary square barrier in the gray area with param- 
eters described in the text. The solid, gray, and dashed lines 
represent the distribution at t — 0, 2, and 5, respectively. 



By use of the above relation and the general product, c 
method, introduced by Gilary et al.\12\, the weightings 
Ck, C n in Eq. lfTT)) can be determined as 



k/n 



J[^(i)]'f(M = 0)& 

/[#/„(*)] V?/»0«0«k 



(21) 



The initial wave packet, ^(x,0), is chosen as a har- 
monic coherent state: 



Oh — 



(22) 



where xq and po denote the position and momentum 
of the center of the packet, while \/2\fb is its mean 
width. For an imaginary square barrier of height Vq = 5, 
the evolution of wave packet is illustrated in Fig|3] with 
(xo,b) = (—10,0.08) and po — 5.23 chosen to be near 
one of the resonance peaks in Figfl] It can be seen that 
the amplitude gets enhanced upon entering the imagi- 
nary barrier due to the resonance. The oscillatory be- 
havior only exists at intermediate time, as depicted by 
the gray line. As the wave packet exits the barrier, the 
amplitudes of both transmitted and reflected waves have 
already been enhanced by 55 and 20 times each. 



III. PT-SYMMETRIC PERIODIC POTENTIAL 

A. Origin of bifurcation 

In the following, we generalize the single imaginary 
potential to a PT-symmetric periodic one. Set the unit 
cell in — d/2 < x < d/2 and denote the eigenfunction by 



$(x) = 4>r,k(x) +A(j> hk (x) 



(23) 



where 4> T} k{x) and <pi^{x) are the right- and left-moving 
scattering states in the unit cell similar to Eq.@. Ac- 
cording to the Bloch's theorem, the wavcfunction in its 



previous cell —3d/2 < x < —d/2 can be represented as: 

^(z) = e- lKd [^ k (x + d) + A^ k (x + d)] (24) 

where K is the Bloch wave number. Matching the bound- 
ary conditions at x = —d/2 gives 



-iKd, 



v,fe(o) = # 



-iKd 1 1 



k ^)=A[e- iK %, k {\)-^, k {~)} 

(25) 



where <p' denotes d<j)/dx. By use of the equality [</fv,fc</>; k ~ 
0, the above equations require the fol- 



<?V,fc<W,fcJlx=-d/2 

lowing condition in order to have a nonzero A: 



U x -U 2 



2[4> r ,h4>'l,k - K^k] \x=d/2 



cos Kd 



(26) 



where 



and the definition of U2 is similar but with the subindiccs 
r and / interchanged. Equation (|26p gives us the disper- 
sion relation for the spectrum. 

Let us write down the form of 4> r j at x = ±d/2 explic- 
itly 



<f>r,k( X ) 
4>l,k{ x ) 



R r e 



-ikx 



Te l 



Te 

„ — ik 



—ikx 



Rie 1 



(27) 



where the amplitudes R 1 T may vary with k. Inserting 
them in Eq. ([26|) gives the relation between the energy 
E = k 2 /2m and K for a general periodic potential: 



T — R r R l 



2T 



l e ikd 



1 



— e~ lKa = cos Kd. 
2T 



(28) 



The fact that we are interested at the branch points 
of PT-symmetric periodic potentials allows us to further 
simplify Eq. (|28]) . First, take the complex conjugate of 
the Schrodingcr equation for </> r 



2^^ + V * {x) 



If we choose the unit cell to also exhibit the PT- 
symmetry, the potential V*(x) will become identical to 
V(—x). Furthermore, E* is the same as E since the 
eigenenergy is real at the branch point. These two prop- 
erties combined with Eq. (f2"5|) tell us that (f>* k (x) must be 
a linear combination of the two scattering states: 



</>*,fc(z) 



k,k(~x) + C<f>r tk (-x). 



(30) 
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This relation gives four equations, each of which corre- 
sponds to matching the right- and left-moving part of 
the wavefunctions on either side of the unit cell: R* = B, 
1 = BRi+CT, = BT+CR r and T* = C. It is straight- 
forward to show that they yield T* = T/(T 2 - R r Ri) 
which simplifies Eq. ([28|) to: 



1 
W\ 



cos(kd + 9) = cos Kd 



(31) 



where 9 is the phase of T. If |T| should ever be larger 
than 1 near the extrema of the LHS in Eq. (|3~lj) . there 
will be no real solution for those K that make | cos Kd\ 
large. In other words, the eigenenergies near the Bril- 
louin zone boundaries of K are sure to be complex. We 
emphasize that the condition that \T\ exceeds unity for 
some k alone is not enough to predict the bifurcation. 
It has to coincide with the occurrence of extreme values 
for cos(kd + 9) /\T\. This is complement to the previ- 
ous conclusion [l0( that perturbation forbids bifurcation. 
What we have shown here is that \T\ < 1 guarantees 
there will be no bifurcation. It is then natural to ask 
whether \T\ < 1 implies the validity of perturbation. As 
far as the PT-symmetric periodic square potential in the 
next subsection is concerned, this statement has been 
checked to hold in Section 11(B). 



B. PT-symmetric periodic square potential 



+ iVo 



a 



-iVo 



FIG. 4: A PT-symmetric periodic potential with its unit cell 
highlighted by the dashed line. Parameters are set at Vo = 5, 
a = 1 and m = 1. 

As a demonstration, we quantify our conclusions in 
the previous subsection for the special case of a PT- 
symmctric periodic square potential, as shown in FigSl 
The transmission coefficient for the unit cell can be cal- 
culated as: 



T = 



i(k + +k~ -2k)a 



1 — r a rb€. 2lk+a + rbr c e 2lk a — r a r c e 2t< -- k++k )° 



(32) 



where k = \j2mE and k^ = yj2vn(E =p iVq) represents 
the momentum in the imaginary barrier/well region. The 
transmission and reflection amplitudes corresponding to 




(a) 
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Bloch wave number [K] 



Bloch wave number [K] 



FIG. 5: (a) The RHS of Eq.(j3TJ, cas{kd + 0)/\T\, is shown in 
the solid line for the potential and parameters defined in Figf3] 
There will be allowed (unshaded) and forbidden (shaded) re- 
gions of k. When cos K d intersects the point highlighted 
by the dashed circle, the corresponding K value illustrates 
one branch point, (b) The real part of the first dispersion 
band. The capital letter R/C indicates that the spectrum is 
real/complex in magnitude, (c) The imaginary part of the 
first dispersion band. 



the three boundaries are denoted by t a = 2k/(k + fc + ), 
t b = 2k+/{k+ + k~), t c = 2k~/(k- + k), r aib = 1 - t a . b 
and r c = t c — 1. 

Inserting Eq. (f3"2"| into Eq. (|3"Tj) and the dispersion rela- 
tion can be determined as: 

(k + — k~) 2 

cos(fc + + k~)a — ; , sin k + a sin k~a = cos 2Ka. 



2k+k- 



(33) 



The band structure is solved numerically and plotted in 
FigEl 

The perturbation method has been demonstrated to 
break down for a single localized complex barrier poten- 
tial in Section 11(B), either with or without an accom- 
panying well. When it is generalized to a PT-symmetric 
periodic potential, we are sure that the perturbation re- 
mains at fault from the analysis within one unit cell. 
In the meantime, Fig[5] shows that the originally real 
eigenenergies evolve into both complex and real spectra. 
So the existence of branch points and the breakdown of 
perturbation are tied together as proved by Ref . (Io| . In 
the mean time, the values of \T\ at these extrema also 
become greater than unity, as is shown in TabHTl 

In passing, we observed a special bond between the 
branch points and the Fano resonances discussed in Sec- 
tion 11(B). Since the momenta k + and k~ here happen 
to be complex conjugates when E is real, the reflection 
amplitude in Eq. (|14p is rendered pure imaginary. The 
condition for constructive interference becomes identical 
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E 


10.39 


30.62 


60.34 


99.87 


\T\ 


1.36 


1.02 


1.00 


1.00 


71 + 1/2 


1.50 


2.50 


3.50 


4.50 



TABLE II: Partial list of energies and absolute magnitude 
of the transmission amplitude at the branch points. Their 
corresponding resonance numbers n from n + 1/2 = 2a/ A = 
a\/2mE /it indeed are close to integers. 



to Eq. p3p . This is supported by the numerical results 
for the branch points presented in TabHTl which indeed 
approach those for the resonances. 



IV. DISCUSSIONS AND CONCLUSIONS 

Since the non-Hcrmitian potential results from our ex- 
amining an open system which is allowed to exchange 
particles with its environment, the conventional continu- 
ity equation is no longer obeyed fill. fiHl [Til. [i~o| . One in- 
teresting debate [l6| is whether it is meaningful to mend 
this loss of unitarity. Several studies 0, Hi [H have 
been devoted to this effort. They tried to construct a 
metric and correspond ing transformed wave function as 
the general framework [l^ US HH, HH, HH to examine the 
subsystem in a quasi- Hermitian analysis. Although the 



unitarity can be restored, one is constrained to discuss 
wavefunctions either coming simultaneously from both 
sides [ill [TtJ or exhibiting different normalization con- 
stants on either side of the potential [l||. To resolve these 
dilemmas, Jones [17] suggested that one simply treats the 
non-Hcrmitian scattering potential as an effective one in 
the standard framework of quantum mechanics. This is 
the view we adopt in this work at studying the localized 
NH potentials as a phenomcnological model of a quantum 
sink or source in an open system. 

In conclusion, although the breakdown of perturbation 
approach has been identified to be intimately linked to 
the existence of complex spectrum for a PT-symmctric 
periodic potential, we offer new insights into their re- 
lationship by studying the less sophisticated case of an 
imaginary barrier potential. Due to the simplicity of this 
potential form, we are able to elucidate the cause of the 
similar breakdown of perturbation, why the perturbation 
can be remedied by the superimposition of a sufficiently 
large real potential, and the condition for a constructive 
interference which differs from that for a real potential. 
Most importantly, we provide a more comprehensive cri- 
teria through Eq. (|3ip on why, when and where the bifur- 
cation will occur. 

We thank Hsiu-Hau Lin for useful comments and ac- 
knowledge the support by NSC in Taiwan under grant 
No. 95-2112-M007-046-MY3. 



[1] C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80, 
5243 (1998). 

[2] C. M. Bender, Contemp. Phys. 46, 277 (2005). 

[3] A. Mostafazadeh, Ann. Phys. (N.Y.) 309, 1 (2004). 

[4] C. M. Bender, Rep. Prog. Phys. 70, 947 (2007). 

[5] C. M. Bender, D. C. Brody, H. F. Jones, and B. K. Meis- 

ter, Phys. Rev. Lett. 98, 040403 (2007). 
[6] A. Mostafazadeh, Phys. Rev. Lett. 99, 130502 (2007). 
[7] P. Dorey, C. Dunning, and R. Tateo, J. Phys. A: Math. 

Gen. 34, 5679 (2001). 
[8] Z. H. Musslimani, K. G. Makris, R. El-Ganainy, and D. 

N. Christodoulides, Phys. Rev. Lett. 100, 030402 (2008). 
[9] K. G. Makris, R. El-Ganainy, and D. N. 

Christodoulides, and Z. H. Musslimani, Phys. Rev. 

Lett. 100, 103904 (2008). 
[10] S. Klaiman, U. Giinther, and N. Moiseyev, Phys. Rev. 

Lett. 101, 080402 (2008). 
[11] F. Cannata, J-P. Dedonder, and A. Ventura, Ann. Phys. 

(N.Y.) 322, 397 (2007). 



[12] I. Gilary, A. Fleischer, and N. Moiseyev, Phys. Rev. A. 

72, 012117 (2005). 
[13] Z. Ahmed, Phys. Lett. A. 324, 152 (2004). 
[14] Z. Ahmed, C. M. Bender, and M. V. Berry, J. Phys. A: 

Math. Gen. 38, L627 (2005). 
[15] M. Znojil, J. Phys. A: Math. Gen. 39, 13325 (2006). 
[16] H. F. Jones, Phys. Rev. D. 76, 125003 (2007). 
[17] H. F. Jones, Phys. Rev. D. 78, 065032 (2008). 
[18] M. Znojil, Phys. Rev. D. 78, 025026 (2008). 
[19] C. M. Bender, D. C. Brody, and H. F. Jones, Phys. Rev. 

Lett. 89, 270401 (2002); ibid. 92, 119902(E) (2004). 
[20] A. Mostafazadeh, J. Math. Phys. (N.Y.) 43, 205 (2002). 
[21] H. F. Jones, J. Phys. A: Math. Gen. 38, 1741 (2005). 
[22] A. Mostafazadeh, J. Phys. A: Math. Gen. 38, 6557 

(2005) ; ibid. 38, 8185(E) (2005). 

[23] H. F. Jones and J. Mateo, Phys. Rev. D 73, 085002 

(2006) . 



